Disappearance of back-to-back high px hadron correlations in central Au+Au 

collisions at ^/snn = 200 GeV 
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Azimuthal correlations for large transverse momentum charged hadrons have been measured over 
a wide pseudo-rapidity range and full azimuth in Au+Au and p+p collisions at ^/snn = 200 GeV. 
The small-angle correlations observed in p+p collisions and at all centralities of Au+Au collisions 
are characteristic of hard-scattering processes already observed in elementary collisions. A strong 
back-to-back correlation exists for p+p and peripheral Au + Au. In contrast, the back-to-back 
correlations are reduced considerably in the most central Au+Au collisions, indicating substantial 
interaction as the hard-scattered partons or their fragmentation products traverse the medium. 

PACS numbers: 25.75 



In collisions of heavy nuclei at high energies, a new 
state of matter consisting of deconfined quarks and glu- 
ons at high density is expected Q . Large transverse mo- 
mentum partons in the high-density system result from 
the initial hard scattering of nucleon constituents. Af- 
ter a hard scattering, the parton fragments to create a 
high energy cluster (jet) of particles. A high momen- 
tum parton traversing a dense colored medium is pre- 
dicted to experience substantial energy loss § U and 
may be absorbed. Measurement of the parton fragmen- 
tation products after hard-scattering processes in nuclear 
collisions may reveal effects due to the interaction of high- 
momentum partons traversing the medium, thereby mea- 
suring the gluon density of the medium 

Hard scattering processes have been established at high 
transverse momentum (j>t) m elementary collisions at 
high energy through the measurement of jets H ^, 
correlated back-to-back jets (di-jets) |l, high px single 
particles, and back-to-back correlations between high pt 
hadrons Q . Jets have been shown to carry the momen- 
tum of the parent parton fill. The jet cross sections and 



high pt single particle spectra are well described over 
a broad range of energies JlT[ in terms of the hadron's 
parton distributions, hard parton scattering treated by 
perturbative QCD, and subsequent fragmentation of the 
parton. High pt jet events have also been studied in 
proton- nucleus interactions |I2| |. In the absence of ef- 
fects of the nuclear medium the rate of hard processes 
should scale linearly with the number of binary nucleon- 
nucleon collisions. Recent results from RHIC, however, 
show a suppression of the single particle inclusive spectra 
of hadrons for px > 2 GeV/c in central Au+Au collisions, 
indicating substantial in- medium interactions Jl3| , [0J |. 

In this Letter, we report measurements of two-hadron 
angular correlations at large transverse momentum for 
p+p and Au+Au collisions at ^/JWn — 200 GeV. These 
correlation measurements provide the most direct ev- 
idence for production of jets in high energy nucleus- 
nucleus collisions, and allow for the first time measure- 
ments, inaccessible in inclusive spectra, of the fate of 
back-to-back jets in the dense medium as a function of 
the size of the overlapping system. The results reveal 
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significant interaction of hard-scattered partons (or their 
fragmentation products) in the medium, with a strong 
dependence on the geometry and distance of traversal. 

The measurements were made using the STAR detec- 
tor lla] at the Relativistic Heavy-Ion Collider (RHIC) 
at Brookhaven National Laboratory. The STAR detec- 
tor is a large acceptance magnetic spectrometer, with a 
large volume Time Projection Chamber (TPC) inside a 
0.5 Tesla solenoidal magnet. The TPC measures the tra- 
jectories of charged particles and determines the particle 
momenta. The TPC has full azimuthal coverage over 
a pseudo-rapidity range \rj\ < 1.5. STAR has excellent 
position and momentum resolution, and, due to its ver- 
texing capabilities, is able to identify many sources of 
secondary particles. The p+p analysis uses «10 million 
minimum bias p+p events triggered on the coincidence 
of signals from scintillator annuli spanning the pseudo- 
rapidity interval 3.5 < |?7| < 5.0. The Au+Au analysis 
uses f»1.7 million minimum bias Au+Au events and «1.5 
million top 10% central Au+Au events. 

Partons fragment into jets of hadrons in a cone around 
the direction of the original hard-scattered parton. The 
leading hadron in the jet tends to be most closely aligned 
with the original parton direction ]l6| ] . The large multi- 
plicities in Au+Au collisions make full jet reconstruction 
impractical. Thus, we utilize two-particle azimuthal cor- 
relations of high pt charged hadrons (l7) to identify jets 
on a statistical basis, with known sources of background 
correlations subtracted. 

Events with at least one large transverse momentum 
hadron (4 < p% ia < 6 or 3 < p% ig < 4 GeV/c), de- 
fined to be a trigger particle, are used in this analysis. 
For each of the trigger particles in the event, we incre- 
ment the number N(A<f>, Arf) of associated tracks with 
2 GeV/c < pt < P% %9 as a function of their azimuthal 
(A</>) and pseudo-rapidity (A?y) separations from the trig- 
ger particle. We then construct an overall azimuthal pair 
distribution per trigger particle, 



• opposite sign a same sign 



D{A4>) = -r^ / dArjN (Acf>, Arj) , 

J * trigger € J 



(i) 



where N tr igg e r is the observed number of tracks satisfying 
the trigger requirement. The efficiency e for finding the 
associated particle is evaluated by embedding simulated 
tracks in real data. In order to have a high and con- 
stant tracking efficiency, the tracks are required to have 
I77I < 0.7, which translates to a relative pseudo-rapidity 
acceptance of I A?7 1 < 1.4. The single track reconstruction 
efficiency varies from 77% for the most central Au+Au 
collisions to 90% for the most peripheral Au+Au colli- 
sions and p+p collisions. 

Identical analysis procedures are applied to the p+p 
and Au+Au data. Displayed in Figure |l| are the az- 
imuthal distributions for same-sign and opposite-sign 
charged pairs from the a) p+p data and b) minimum 
bias Au+Au data for 4 < p% 19 < 6 GeV/c. The data 
are integrated over the relative pseudo-rapidity range 
< |A?7| < 1.4. Clear correlation peaks are observed 
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FIG. 1: Azimuthal distributions of same-sign and opposite- 
sign pairs for a) p+p, b) minimum bias Au+Au, and c) 
background-subtracted central Au+Au collisions. All corre- 
lation functions require a trigger particle with 4 < p^T 13 < 6 
GeV/c and associated particles with 2 GeV/c < pr < V^p 3 ■ 
The curves are one- or two- Gaussian fits. 



near A<f> ~ and A<p ~ it in the data. The opposite- 
sign correlations at small relative azimuth are larger than 
those of the same-sign particle pairs, while the sign has 
a negligible effect on the back-to-back correlations. 

To isolate the jet-like correlations (localized in A</>, 
A77) in central Au+Au collisions, the azimuthal distri- 
butions are measured for two regions of relative pseudo- 
rapidity, |A?7| < 0.5 and 0.5 < \Arj\ < 1.4 0. The 
difference between the small and large relative pseudo- 
rapidity azimuthal distributions is displayed in Figure |l|c 
along with single Gaussian fits. Near A<fi = 0, the az- 
imuthal distributions from Au+Au and p+p have simi- 
lar shapes. For the opposite-sign azimuthal distributions, 
the Gaussian widths are 0.17±0.01(siai.)±0.03(s?/s.) ra- 
dians for p+p data, and 0.20 ± 0.02(stat.) ± 0.03(sys.) 
radians for the central Au+Au data. For the same- 
sign azimuthal distributions, the Gaussian widths are 
0.16 ± 0.02(stat.) ± 0.03(sys.) radians for p+p data, and 
0.15 ± 0.03(stat.) ± 0.04(sys.) radians for the central 
Au+Au data. The systematic errors reflect the spread 
of values found for different choices of the A</> bin width. 
Within the errors, there are no significant differences be- 
tween the small-angle correlation widths for p+p and 
central Au+Au collisions. 

The ratios of the opposite-sign to same-sign peak ar- 
eas are 2.7 ± 0.9(stat.) ± 0.2(sys.) for p+p and 2.5 ± 
0.6(stat.) ± 0.2(sys.) for central Au+Au collisions. In jet 
fragmentation, there are dynamical charge correlations 
between the leading and next-to-leading charged hadrons 
]l8| that originate from the formation of qq pairs along 
a string between two partons. This results in a preferen- 
tial ordering into oppositely-charged adjacent particles 
along a string during fragmentation. The Hijing event 
generator, which utilizes the Lund string fragmentation 
scheme |ll| incorporating these concepts, predicts a ratio 
of 2.6 + 0.7 for the opposite-sign to same-sign correlation 
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strengths. The agreement of this ratio with those mea- 
sured in the central Au+Au and p+p suggests that the 
same jet production mechanism is responsible for a ma- 
jority of the charged hadrons with pr > 4 GeV/c in p+p 
and central Au+Au collisions. 

The decay of resonances would also lead to small-angle 
azimuthal correlations, but a resonance decay origin is 
unlikely due to the observed correlation of particles with 
the same charge sign, the similarity in the measured 
small-angle azimuthal correlation widths in the Au+Au 
and p+p interactions, and the strong back-to-back corre- 
lations of large pr particles seen for p+p collisions in Fig. 
|l|a. The latter correlations, indicative of di-jet events H, 
are removed from the central Au+Au sample by the sub- 
traction in Fig. 0c. A quantitative analysis of back-to- 
back jet survival in Au+Au requires the more detailed 
treatment of background correlations described below. 

In addition to correlations due to jets, the two-particle 
azimuthal distributions in Au+Au exhibit a structure at- 
tributable to an anisotropy of single particle production 
relative to the reaction plane. Previous measurements 
p7[ indicate that, even at large transverse momentum, 
the particle distributions contain an anisotropy due to 
elliptic flow that can be characterized by dN/d(<fi—§ r ) oc 
1 + 2V2 cos(2(0— 3?r)), where <f> r is the reaction plane an- 
gle determined event by event and i>2 is the elliptic flow 
parameter. This leads to a two particle azimuthal distri- 
bution of the form, dN/dA<j) = B(l+2«| cos(2A</>)). The 
elliptic flow component of the two-particle azimuthal dis- 
tribution is measured using several methods |jl7| . In this 
paper, V2 is determined using a reaction-plane method. 

A simple reference model can be constructed for the 
two-particle azimuthal distributions of high px particles 
in Au+Au collisions. A number of independent hard 
scatterings (each similar to one measured in a triggered 
p+p event) included in an event with correlations due to 
elliptic flow can be represented by the azimuthal distri- 
bution, 



D 1: 



odd 



D pp + B(l + 2v\ cos(2A<?(>)). 



(2) 



The elliptic flow parameter (v^) is measured indepen- 
dently in the same set of events, and is taken to be con- 
stant for pt > 2 GeV/c |l7| . The parameter B is then 
determined by fitting the observed £) AuAu in the region 
0.75 < \A<f>\ < 2.24 radians, which is largely free of jet 
contributions in the p+p data. 

In Figure]^, the azimuthal distributions for < \Arj\ < 
1.4 in Au+Au collisions at various centralities arc com- 
pared to Equation ^ using the measured p+p data. The 
centrality selection is constructed by subdividing the 
Au+Au minimum bias data sample into subsamples with 
different charged particle multiplicities within |7j| < 0.5. 
The parameters v-i and B are determined independently 
for each centrality bin, and are listed in Table |. For 
all centralities, the azimuthal correlation near A(j) = 
is well described by Equation ^. This indicates that the 
same mechanism (hard parton scattering and fragmenta- 
tion) is responsible for high transverse momentum parti- 
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FIG. 2: Azimuthal distributions (0 < At;| < 1.4, 4 < p™ 9 < 
6 GeV/c) for Au+Au collisions (solid circles) compared to the 
expected distributions £) modcl from Equation^ (open circles). 
Also shown is the elliptic flow contribution for each centrality 
(solid curve). 



Centrality (%) 


N part 


V2 


B 


60-80 


20+6 


0.24+0.04 


0.065+ 0.003 


40-60 


61+10 


0.22+0.01 


0.231+ 0.003 


30-40 


114+13 


0.21+0.01 


0.420+ 0.005 


20-30 


165+13 


0.19+0.01 


0.633+ 0.005 


10-20 


232+11 


0.15+0.01 


0.931+ 0.006 


5-10 


298+10 


0.10+0.01 


1.187+ 0.008 


0-5 


352+7 


0.07+0.01 


1.442+ 0.003 



TABLE I: Centrality, number of participants, V2 (2 < pr < 6 
GeV/c), and normalization constant B. The errors on V2 
and B are statistical only, while the errors on the number of 
participants are systematic |14|. 



cle production in p+p and Au+Au collisions. However, 
the back-to-back correlations are suppressed in Au+Au 
collisions compared to the expectation from Equation ||, 
and the suppression is greater for more central collisions. 
The most central collisions show no indication of any 
back-to-back correlations beyond that expected from el- 
liptic flow. 

The ratio of the measured Au+Au correlation excess 
relative to the p+p correlation is: 



iAA (A0i,A<£ 2 ) 



-B(l+2vl cos(2A0))] 



(3) 



The ratio can be plotted as a function of the number of 
participating nucleons {N part ), deduced from the central- 
ity bins as described in reference Jl4| . Iaa is measured 
for both the small-angle (|A</>| < 0.75 radians) and back- 
to-back (|A</>| > 2.24 radians) regions. The ratio should 
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FIG. 3: Ratio of Au+Au and p+p (Equation for small- 
angle (squares, |A0| < 0.75 radians) and back-to-back (cir- 
cles, \A<f>\ > 2.24 radians) azimuthal regions versus num- 
ber of participating nucleons for trigger particle intervals 
4 < pp 9 < 6 GeV/c (solid) and 3 < V % %a < 4 GeV/c (hol- 
low). The horizontal bars indicate the dominant systematic 
error (highly correlated among points) due to the uncertainty 
in V2- 



be unity if the hard-scattering component of Au+Au col- 
lisions is simply a superposition of p+p collisions unaf- 
fected by the nuclear medium. These ratios are given 
in Figure |^ for the trigger particle momentum ranges 
indicated. The asymmetric systematic errors are domi- 
nated by the +5/-20% systematic uncertainty on vi due 
to the potential non-flow contributions [POf as well as 
other sources of systematic uncertainty EW T 

For the most peripheral bin (smallest N part ), both 
the small-angle and back-to-back correlation strengths 
are suppressed compared to the expectation from Equa- 
tion g This may be an indication of initial state nu- 
clear effects such as shadowing of parton distributions 
or scattering by multiple nucleons, or may be indica- 
tive of energy loss in a dilute medium pM. As N part 
increases, the small-angle correlation strength increases, 
with a more pronounced increase for the trigger parti- 
cles with lower px threshold. If there were a large non- 
jet contribution to particle production (i.e. collective 
transverse flow) at the trigger threshold and above, it 
would dilute the jet related correlation signal and this 
ratio would be reduced. The back-to-back correlation 
strength, above background from elliptic flow, decreases 
with increasing N part and is consistent with zero for the 
most central collisions. In the extreme case, if there 
were no elliptic flow for the 0-5% most central collisions, 
I AA (2.24:,n) = 0.4 + 0.1 for 4 < p™ 9 < 6 GeV/c, com- 
pared to Iaa{2-24: 1 7t) = 0.1 ± 0.1 using the measured 
elliptic flow value. Therefore, an overestimation of the 
elliptic flow cannot explain the observed suppression of 
back-to-back correlations. 

Analyses of fixed-target experiments have sug- 
gested that the shape of the back-to-back dihadron az- 
imuthal distribution is sensitive to the intrinsic parton 
transverse momentum kx in the initial state. In proton- 



nucleus and nucleus- nucleus collisions, additional initial- 
state transverse momentum can be generated by multi- 
ple nuclcon-nucleon interactions preceding a hard scat- 
tering (|^, [24], [29]. To investigate whether this nuclear 
hr can account for the observed deficit of back-to-back 
azimuthal correlations in central Au+Au collisions, we 
have carried out Pythia ]l9| simulations varying the kx 
parameter. A rather extreme change from the nominal 
value of (7 = 1 GeV/c to 4 GeV/c introduced only a 
small effect, reducing the predicted 7^(2.24,71-) by less 
than 20%. Experimental study of initial state effects on 
the azimuthal correlations requires the measurement of 
p(d)+Au collisions at RHIC energies. 

In addition to the present data, two other striking 
effects have been observed at high px in nuclear colli- 
sions at RHIC: strong suppression of the inclusive hadron 
yield in central collisions M, and large elliptic flow 
which saturates at px > 3 GeV/c |i7| |. These phenom- 
ena are all consistent with a picture in which observed 
hadrons at px > 3 — 4 GeV/c are fragments of hard scat- 
tered partons, and partons or their fragments are strongly 
scattered or absorbed in the nuclear medium. The ob- 
served hadrons therefore result preferentially from hard- 
scattered partons generated on the periphery of the reac- 
tion zone and heading outwards |2q |. In this picture the 
inclusive yield will be suppressed relative to the binary 
scaling expectation, and the strong position-momentum 
correlation required to explain the large elliptic flow |27j] 
emerges naturally. The properties of small- angle hadron 
correlations will have weak dependence on the size of the 
colliding system, whereas the back-to-back correlations 
will exhibit strong suppression for a large system relative 
to a small one, both as observed. 

In summary, STAR has measured azimuthal correla- 
tions for high px charged particles over a large relative 
pseudo-rapidity range with full azimuthal angle coverage. 
Comparison of the opposite-sign and same-sign correla- 
tion strengths indicates that hard scattering and frag- 
mentation is the predominant source of charged hadrons 
with px > 4 GeV/c in central Au + Au collisions. The 
azimuthal correlations in Au+Au collisions have been 
treated as the superposition of independently determined 
elliptic flow and individual hard parton scattering con- 
tributions, the latter measured in the STAR p+p data. 
The most striking feature of the hard-scattering com- 
ponent is an increasing suppression of back-to-back rela- 
tive to small-angle correlations with increasing centrality. 
These observations appear consistent with large energy 
loss in a system that is opaque to the propagation of high- 
momentum partons or their fragmentation products. 
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